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■ 1. Introduction 

o 



Abundances of chemical elements in galaxies are commonly measured using the emission lines 
emitted by astrophysical nebulae (e.g., see Osterbrock 1989). Recombination lines of hydrogen 



2 ■ (e.g., the Balmer series) or helium and forbidden lines of singly and doubly ionized carbon, nitro- 

gen, oxygen, neon, and sulfur are among the strongest observable lines in the visible and ultraviolet 
spectrum. The relative strengths of each emission line, combined with some knowledge of the 



temperature and density of the nebulae, provide information about the relative and absolute abun- 
dances of each ion. Chemical analyses of individual nebulae have been used to test nucleosynthesis 
models of stars, the chemical evolutionary history of galaxies, and nucleosynthesis in the big bang 
(review by Aller 1990; Pagel 1998). Even integrated spectra of entire galaxies may be used to 
estimate the overall degree of chemical enrichment of a galaxy (Kobulnicky, Kennicutt, & Pizagno 
1998). 

Observational techniques sometimes permit the relative fluxes of emission lines from galaxies to 
be measured with a high degree of precision. However, in the current generation of wide-field galaxy 
surveys on multi-object spectrographs, flux calibration is frequently problematic due to unfavorable 
observing conditions or instrumental effects such as a variation in system response over the field 
of view. In addition, the traditional line ratio for reddening determination is Ha/ HP, but for 
redshifts greater than z ~ 0.3 the Ha line is not available, and these galaxies rarely have sufficient 
signal-to-noise in their spectra to attempt a reddening derivation from the higher-order Balmer 
lines. For such surveys, it is still desirable to extract chemical information from the data, where 
possible. 
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In this paper we explore the possibility of estimating mean gas-phase oxygen abundances for 
galaxies based the equivalent widths (rather than fluxes) of strong Balmer HP recombination line 
and forbidden [0 II] and [0 III] emission lines. This method essentially uses the underlying stellar 
continuum as a crude flux calibrator, and has the advantage that any reddening affecting the stars 
and gas equally will be automatically removed. While we would not expect the equivalent-width 
method to be as accurate as the traditional method, particularly for individual objects, we find 
that it works sufficiently well to enable us to measure approximate oxygen abundances in relatively 
faint, intermediate-redshift galaxies. 

Our approach is to develop the method from basic principles, and to test and empirically 
calibrate it using three spectroscopic galaxy datasets, and in the end compare the oxygen abundance 
indicator R23 to the quantity EWR23 derived from traditional line fluxes and flux ratios to results 
using only the equivalent widths and equivalent width ratios. The thesis of this paper is that even 
when the dereddened line flux ratio R23 is not available, a corresponding ratio of equivalent widths 
can still provide an estimate of the oxygen abundance. 

This new method is applied to the study of abundances in a set of intermediate-redshift galaxies 
observed as part of the Deep Extragalactic Evolutionary Probe (DEEP) spectroscopic survey in 
the Groth Strip (DGSS) presented in a companion paper, Kobulnicky et al. (2003; Ke03). 

2. Gas-phase Oxygen Abundance Measurements from Emission Lines 

Standard techniques for measuring chemical abundances in astrophysical nebulae from the 
fluxes and flux ratios of nebular emission lines are reviewed in Osterbrock (1989). The most direct 
and reliable techniques involve measuring a suite of temperature-sensitive and density-sensitive 
line ratios to determine the physical conditions of the plasma. In the absence of high signal- 
to-noise data measuring temperature-sensitive line ratios (e.g., [O III] A4363/[0 III] A5007), the 
default diagnostic for measuring the oxygen abundance of ionized nebulae has become the ratio of 
strong oxygen emission lines R23 = (I3727 + ^4959 + ^5007) (Pagel et al. 1979). Subsequently, 
many authors have developed formulations relating this strong line ratio to the gas-phase oxygen 
abundance (e.g., Edmunds & Pagel 1984; McCall, Rybski & Shields 1985; Dopita & Evans 1986; 
McGaugh 1991). Most modern calibrations relating R23 to oxygen abundance include a measure of 
the ionization parameter, such as O32 = (/4959 + ^5007)/-^3727 = I\o ni]l 'I[o 11] as a second parameter 
(e.g., McGaugh 1991; Pilyugin 2001). 

While the R23 method was developed primarily for individual HII regions, Kobulnicky, Ken- 
nicutt, & Pizagno (1998; KKP) have shown that this approach is valid even when applied to the 
integrated spectra of galaxies as a means of estimating an average O/H abundance for the galaxy as 
a whole. Even when the signal to noise ratio of these lines is as low as 8:1 or when a spectroscopic 
observation encompasses a range of metallicities and ionization conditions within a galaxy, this ra- 
tio provides a rough, but reliable, estimate of the mean gas-phase oxygen abundance when used in 
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conjunction with the appropriate calibration relating R23 to O/H. At this signal-to-noise ratio, the 
associated uncertainty on the oxygen abundance due to line measurement error is typically com- 
parable to the uncertainty due to the calibration of the -R23 vs. O/H relationship: ~ 0.15 dex. In 
this paper we show that the uncertainty introduced by substituting emission line equivalent widths 
is less than these other sources of uncertainty, and thereby establish a new method for measuring 
interstellar medium oxygen abundances. 



3. Justification of the Equivalent Width Method 

The original .R23 method involves line ratios, so we consider the effects of using equivalent width 
ratios in place of flux ratios. For generality, we include the effects of reddening in the discussion, 
assuming the traditional "screen" model. The equivalent width 1 can be written as 

Wx = (1) 

where F\ represents the line flux and Fcx is the underlying continuum flux. The fluxes do not need 
to be calibrated, since the calibration affects both quantities equally. The dereddened, calibrated 
flux value for a line, I\, or continuum, F®, is given by 

/ A = i^ = J F A xlO c(1+/(A)) , (2) 

where c is the logarithmic attenuation at Hf3 and /(A) is a function describing the reddening curve 
(see Seaton 1979). Combining these leads to a general line ratio formulation of 

Al = Wxi Fcxi 1QC( / ( A1)-/(A2)) = f3l 
I\2 Wx2F C X2 W X 2 ' ^ ' 

where we group all the continuum and reddening terms into the factor a. Note that the continuum 
will be attenuated by an amount characterized by c*, and reddened by (presumably) the same 
reddening law. Thus, 

a _ F CXi 10 (c-c*)(f(M)-f(X2))_ ^ 

F CX2 

The factor a contains two unknowns, the ratio of the non-reddened continuum fluxes, which 
depends on the underlying stellar population, and the difference in the attenuation of the emitting 
gas and the continuum light. Note that values of c can range from zero to quite large for individual 
HII regions, although large values are severely deweighted in the average over an entire galaxy. 
Typical derived values are in the range to 1, with a median value around 0.3 (McCall et al. 1985; 
Olofsson 1995). Values of c* can be estimated from the relation c* = 1A7Eb-v from Seaton. For 



1 We adopt the sign convention that emission line equivalent widths are positive. 
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Eb-v m a typical range of 0-0.2, c* ~ 0-0.3. In physically realistic situations, we would expect 
some correlation of c and c*, and also that (c — c*) > 0. 

Let us now consider the specific line ratios of interest for the i?23 method. The first is 
I\0 in]/ Ihp- For the purposes of reddening correction, we adopt the wavelength of the stronger 
[O III] A5007 line. The value of a in this case is 



7?o 

C4861 



a/33 = ^po7i ( c - c *)(- - 034 ). (5) 



Due to the proximity of these emission lines in wavelength, we expect the flux ratios of the 
underlying stellar population to be very close to unity. Similarly, for realistic < (c — c*) < 1, the 
reddening factor will be between unity and 0.92. Thus, for this line ratio, a ~ 1, and we will ignore 
it in further discussion. To a very good approximation, 

ho in] „ w \oni] _ (g) 
Ih(} Wh/3 

The second ratio, I[oii]/Ih/3, is more problematic; for these lines, 

a2/3 = ^372Il0(^)(O-255). ( 7 ) 
^C4861 

In principle, we could estimate from a detailed spectral analysis of the flux-calibrated, integrated 
stellar spectrum and Balmer line ratios, but for faint, galaxies with redshifts z > 0.3 we will have 
neither the signal-to-noise nor access to the Ha line in order to do this. In practice, we need to 
adopt an appropriate average value of a^p determined by other means. 

Galaxy light tends to be dominated by A main sequence and G and K giants (e.g., Morgan 
& Mayall 1957; Pritchet 1977; Kobulnicky & Gebhardt 2000), and assuming an underlying stellar 
spectrum composed of a linear combination of these two spectral types implies that the dereddened 
flux ratio of the continuum ranges from ~1 for late-B stars to ~0.4 for mid-G giants to ~0.2 
for early-K giants (although galaxies whose A3727 flux is dominated by K stars is unlikely to 
have any line emission from star formation). Examining the spectra in Kennicutt (1992b) shows 
the A3727-to-A4861 ratio of the continuum fluxes in galaxies with obvious emission lines ranging 
from ~0. 4-1.0; if these were dereddened the range would shift upwards and might possibly narrow 
somewhat. 

On the other hand, the reddening correction ranges from 1 to 1.8 for realistic values < (c — 
c*) < 1, with a typical likely difference (c — c*) ~ 0.3 giving a reddening factor of ~1.2. Combining 
the reddening and stellar population factors leads to an expected value of 02/3 ~ 0.84 ± 0.3. This 
expected average value for 02/3 differs from unity by less than 0.1 dex. 

We find the R23 measure can now be expressed 

^3 = log /[Q/J1 + /[om] = log ^3om±3om. (8) 

Ih/3 W H /3 
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Not surprisingly, the R23 measured using equivalent widths and an adopted average value for a 2/ 3 
will be most in error when [O II] dominates [O III]. 

The final line ratio of interest is the ionization parameter, and it is easy to see that 

~ W l° ni ^ . (9) 

I[OII] «23^[0/7] 

The value 023 — 02/3 to a very good approximation, so at this point we will drop the subscripts 
and simply refer to "a" for either of these quantities. Comparing the ratios above is probably 
the easiest means to empirically estimate an average value for a. While it would be equally valid 
to estimate a from the [O II] and H/3 ratios, the H/3 equivalent width can be contaminated by 
underlying stellar absorption. 



4. Data Selection and Analysis 

In order to empirically measure an average value of a and to confirm the validity of using 
equivalent width line ratios in estimating R23, we compiled three sets of emission- line spectra of 
nearby galaxies. The first set consists of 16 objects with spatially-integrated spectra from the 
55-object spectroscopic galaxy atlas of Kennicutt (1992a,b; K92) plus six additional emission- 
line objects from KKP. We refer to this sample as the K92+ sample. The K92+ spectra are 
produced by drifting a longslit across each galaxy and have spectral resolutions of 5-7 A (K92) and 
3 A (KKP). The 16 K92 galaxies in our subsample are the strongest emission-line objects where 
global metallicity measurements are possible. The full set of K92 galaxies represents a range of 
morphological types from Sa to Im, but it includes only the bright galaxies of each type. The 
KKP galaxies are all under luminous, dwarf emission-line galaxies. For a second local sample, 
we examined the 198-galaxy Nearby Field Galaxy Survey (Jansen et al. 2000a,b; NFGS), which 
is selected from the CfA redshift catalog (Huchra et al. 1983). From this spatially- integrated 
spectroscopic survey, we culled 98 objects with measurable Hf3, [O III] and [O II] emission lines. 
These spectra have a resolution of 6 A and include a larger range of luminosity (—14 < Mb < —22) 
and surface brightnesses than K92+ while spanning morphological types. The K92+ objects have a 
higher fraction of star-forming galaxies (objects with strong emission lines) compared to the NFGS 
sample. As a third local sample, we used emission-line selected galaxies in the Kitt Peak National 
Observatory Spectroscopic Survey (KISS; Salzer et al. 2000). KISS is a large-area objective prism 
survey of local (z < 0.09) galaxies selected by strong Ha emission, and thus, preferentially contains 
objects with active star formation. While high-quality slit spectroscopy has been obtained for 
~519 KISS galaxies (Melbourne & Salzer 2002), 396 galaxies lack [O II] measurements, leaving 123 
galaxies. These remaining spectra have 5-8 A resolution, and are analyzed in Salzer et al. (2003). 

In summary, the 22 K92+, 98 NFGS, and 123 KISS galaxies were selected from their larger 
parent samples because of strong emission lines suitable for nebular metallicity measurement. Fol- 
lowing KKP, we chose only galaxies with detectable [O II]A3727, [O III]A4959, [O III]A5007, and 
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H/3 emission lines. Only galaxies where all four emission lines were measured with a S/N of 8:1 or 
greater were retained. This selection criterion preferentially includes galaxies with high equivalent 
width lines, but it also includes galaxies with low equivalent width lines where the continuum is 
smooth and well- measured. 

For each galaxy in the K92+ sample we measured the emission-line fluxes and equivalent 
widths manually using Gaussian fits. For the NFGS and KISS surveys, we adopted the published 
emission-line fluxes and equivalent widths. Dereddened emission-line fluxes were computed for all 
three samples by comparing the observed Faa/Fiip ratios to theoretical ratios. 2 The line fluxes are 
dereddened using the law of Seaton (1979) as parameterized by Howarth (1983) and as described in 
Kobulnicky & Skillman (1996). We did not correct the Balmer emission lines for underlying stellar 
absorption. The effects of Balmer absorption by the stellar population are discussed in Section 5.3 

While the K92+KKP sample consists entirely of starforming galaxies, several low-level AGN 
are known to exist in the NFGS. These four objects are not included in the our subsample. The 
123 KISS galaxies included here do not contain any AGN, as they were selected to be conventional 
starforming galaxies based on analysis high-quality spectroscopic observations (Salzer 2003). In 
any case, the presence of AGN among the samples would not have a significant bearing on the 
results of this paper since we are interested in comparing observable properties of emission lines 
rather than deriving physical quantities such as density or metallicity which are sensitive to the 
nature of the ionizing source. 

5. Analysis of Emission Line Quantities 

5.1. [O II] and [O III] Equivalent Widths versus Fluxes 

Figure 1 compares the oxygen and hydrogen emission-line flux ratios to equivalent width ratios 
as a function of emission line strength and B — V color for the K92+ and NFGS samples. Solid 
symbols denote the Nearby Field Galaxy Sample while crosses denote the K92+ galaxies. The upper 
left panel compares the ratio of dereddened [O II] to H(3 fluxes, I^q u]/Ih/3i versus the ratio of [O II] 
to H(5 equivalent widths, W\q ii\/Whp- A solid line marks the 1-to-l correspondence. There is 
a good correlation between the two quantities, indicating that strong-line equivalent widths are 
a good surrogate for dereddened line fluxes. The RMS deviation from the 1-to-l correspondence 
is &{log[W[o ii]/Wh/3]) = 0.11 dex for the combined K92+NFGS samples. Similarly, the panel 
at top right shows the I\o iii]/Ih/3 versus Wiq iih/Wh/3 ratios and indicates that these equivalent 
width ratios are an excellent substitute for dereddened line fluxes, as expected. The RMS deviation 
from the 1-to-l correspondence is cr(log[W[o iiij/Wup]) = 0.05 dex for the combined K92+NFGS 
samples. 



2 lHa/lH/3 = 2.75-2.86 for wide temperature range, e.g., (Hummer & Storey 1987). Here we assumed fixed electron 
temperature of 12,000 K so that Ihcx/Ihp = 2.85. 
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The lower panels of Figure 1 show residuals from the 1-to-l line as functions of Wup, W[qii], 
the ratio Wiq ni]/W<o jn and galaxy color. The K92+ and NFGS galaxies have very small residuals 
in the right column which compares [O III] to Hf3 ratios. Although small, the NFGS residuals are 
slightly systematic, with an offset of 0.07 dex. The excellent correspondence of [O III] fluxes 
to equivalent widths may be easily understood since the [O III] AA4959,5007 lines are close in 
wavelength to H/3 so that neither changes in the underlying galaxy continuum light nor relative 
extinction will alter this ratio. The same data is shown for the [O II] to HP ratios in the left 
column. Here, the K92+ and NFGS galaxies have small, slightly-systematic residuals, particularly 
as a function of color. In these cases, a is unlikely to be quite near to unity. There seems to 
be a systematic offset between the NFGS and K92+ sample, possibly caused by a systematic 
overestimate of the reddening correction to the [O II] A3727 flux due to low spectral resolution in the 
K92 atlas. Within the NFGS sample, the largest residuals are seen in the reddest galaxies, consistent 
with expectations that the reddest galaxies will have the smallest value of a (see Equation 9). The 
systematic trend seen in the residuals with respect to Wup is most likely a reflection of a correlation 
between Wh/3 and color, i.e., redder galaxies have relatively less star formation and hence relatively 
weaker Wh/3- These same galaxies with low Wh/3 will also be the most affected by lack of correction 
for underlying stellar absorption. Nevertheless, the bulk of the residuals are near zero, indicating 
that a has a typical value near unity. 

Figure 2 shows the same comparison as Figure 1 but now for the KISS and K92+ samples. 
There appear to be somewhat more scatter in the KISS sample than in NFGS, and here the RMS de- 
viation from the 1-to-l correspondence is <r(log[IF[o ii]/Wh[3]) = 0.15 dex and <r(log[VF[o ih]/Whp\) = 
0.04 dex for the combined K92+KISS samples. 

The lower panels of Figure 2 show residuals from the 1-to-l line as a function of equivalent 
width and galaxy color. The right column shows very small and non-systematic residuals, indi- 
cating excellent correspondence between [O III] equivalent widths and fluxes for all KISS galaxy 
colors and line ratios. However, the left column shows significant dispersion of 0.15 dex between 
log^o^/H 7 "^] and log[/[£ ii\/Ihp\- The residuals are not correlated with galaxy color or equiv- 
alent width, suggesting that measurement errors and/or uncertainties in the reddening correction 
are responsible for the dispersion. The [O II]/H/3 residuals appear slightly offset from the zero line, 
indicating a slightly above unity on the average. 

5.2. The Ionization Parameter [O III]/[0 II] 

We consider next the ratio of the [O II] and [O III] lines, as this was shown to be the most 
direct way to estimate the unknown value a. In addition, the quantity I\o in]/ I\o II] * s a measure 
of the hardness" of the ionizing photons, and is used in the R23 to O/H conversion. Modern 
calibrations relating R23 to oxygen abundance include a measure of the ionization parameter, 
usually 32 = (I 4 959 + ^5007)/^3727 = I[o in] I 1 [O II] as a second parameter (e.g., McGaugh 1991; 
Pilyugin 2001). 
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In Figure 3, we compare the equivalent width and flux ratios for the K92+ galaxies (crosses) 
and NFGS galaxies (solid symbols). For completeness sake, on the left we plot W\q jj n j W\q jn 
against the observed F\o ni]l F[o u] ratio; there is a clear correlation, but the points are scattered 
mostly below the 1-to-l correlation (shown by the line), consistent with the lack of reddening 
correction. Residuals to the 1-to-l line are shown in the lower panels, plotted against other galaxy 
parameters The strongest correlation is seen against color, which is completely expected: the ratio 
(W[oiii]/W[ OII ])/(F[ OIII ]/F[ OII ]) is identical to the ratio of the continua, ^3727/^4959,5007, 
which will be strongly correlated with B — V. Most of the other systematic trends in the residuals 
can be explained as correlations of the other parameters with B — V. 

The right column shows a comparison of equivalent width ratios to the dereddened line ratios. 
This is of considerable interest, since the ratio (Wiq iin/W\o 11])/ \I\o in\l 'I\o 11]) 1S a direct measure 
of a (see Equation 9). Here there is a much tighter correlation, although there are still some 
systematic trends in the residuals, again most strongly connected to color. We note, however, 
W[o iii]/W[o ii] is actually very close to I[o in]/ I]o ii]-> which confirms that a is near unity on the 
average. We also note that there appears to be a systematic offset between the NFGS and K92+ 
galaxies, perhaps related to the offset seen in the NFGS [O III] data (Figure 1, middle column). 
While we might try to define a color-dependent term for a, it is unclear how much of the systematic 
trend is due to the offset between samples. Without any color-dependent term, the scatter in a 
about unity is only a = 0.12 dex. 

Figure 4 shows the same comparison for the KISS and K92+ samples. Here, the tight relation- 
ship between B — V and the residuals to (W[q ni]/W[o ii])/(F[o 111]/ F[q 11]) is surprisingly lacking. 
Since we are effectively comparing i 7 cA3727/-FcA4959,5007> measured spectroscopically, with the B — V 
color determined from images, this suggests considerable scatter in the continuum measurements. 
Since the KISS data span a redshift range up to < z < 0.1, it is possible that some of the scatter 
is due to [O III] contamination of the V broadband flux, but this would tend to reduce the B — V 
colors of those objects with high Wy ni] at higher z. When we correct for this in a reasonable 
fashion, we do not see a significant reduction in the scatter, so we must assume the scatter due 
to some other reason. The cause is almost certainly that the KISS spectroscopy is not integrated 
across the galaxy, so the spectra tend to be dominated by the nuclear regions, whereas the B — V 
colors are from the integrated galaxy light. 

Looking at the dereddened line ratios in these data, we do not find any trend in a with color, 
although such a trend might be obscured by the scatter. Nevertheless, the scatter about unity is 
only a = 0.14 dex, and again a value of a near unity is indicated. 

In summary, the local data empirically confirm that adopting a = 1 is a reasonable approxi- 
mation, with errors around 35%, and therefore the ratio W]p in] /W]q 11] is a reasonable substitute 
for the ionization parameter, I[o ni]/I[o II]- 
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5.3. EWR 23 versus R23 

The previous sections suggest that the equivalent width ratios can be simply substituted into 
the i?23 formula, and here we provide a direct test of that. We define the measure 

z^w o - W[0 //]A3727 + W[o 777]A4959 + W[o 7//JA5007 n x 
h,W K23 = 777 • UUJ 

which should be approximately equivalent to R23. For completeness, we will compare EWR23 
against both the reddening-corrected R23 and the analogous measure formed with non-corrected 
line fluxes, i2| 3 . Figure 5 (upper left panel) shows the comparison between R^ and EWR23 f° r 
K92+ and NFGS galaxies constructed from the raw fluxes and equivalent widths. A solid line 
illustrates a 1-to-l correspondence. The lower rows in the left column of Figure 5 show residuals 
from the 1-to-l correspondence as a function of Wh/3, W\oii], W\q iii]/W\o ii] and galaxy B — V 
color. The correlation between R^ 3 and EWR23 is strong but has considerable scatter. Formally, 
the RMS dispersion from the 1-to-l relation is a(log[i?23]) = 0-12 dex for the combined K92+NFGS 
samples. 

The correlation is strongest for objects in the K92+ sample and for objects with large values 
of i?23- Deviations from 1-to-l are greatest for objects in the NFGS sample which have low Wnp, 
high W[q ni]/W[o 11] ratios, and red colors. These systematic residuals may be understood as 
a consequence of the lack of corrections for extinction. [O II] A3727 is significantly affected by 
extinction compared to the H/3 and [O III] lines. The measured [O II] flux is a lower limit to 
the true unextincted intensity whereas the measured [O II] equivalent width should be unaffected 
by extinction provided that the extinction toward the gas and stars are similar (but see Calzetti, 
Kinney & Storchi-Bergmann 1994 for evidence that this assumption is sometimes invalid). The 
strong systematic residuals with B—V color seen in the lower left panel is a most likely a consequence 
uncorrected extinction, since redder galaxies are often those with greater extinction. 

In the right column of Figure 5 we show a similar comparison of EWR23 with R23 , where R23 
has been corrected for reddening using the theoretical Balmer decrement. Here the correlation is 
much stronger. The strong correlation between galaxy color and residuals seen in left column is now 
mostly gone, suggesting that the [O II] line fluxes have been successfully corrected for reddening. 
The RMS dispersion from the 1-to-l relation is <j(log[i?23]) = 0.07 dex for the combined K92+NFGS 
samples. The ratio of equivalent widths is a good substitute for the reddening-corrected R23 ratio. 
Use of equivalent widths will be superior to line ratios if the reddening corrections are not known, 
as in the case of galaxies for which the Hq/H/3 ratio is not available (typically true for redshifts 
z > 0.3). The residuals in the lower panels of column 2 are mostly symmetric about zero, with the 
largest scatter again occurring for objects with very low values of Whs- Some of the systematic 
residuals are probably also caused by varying continuum shapes, especially among the NFGS, which 
affect the equivalent widths of the [O II] lines in a systematic manner which is related to galaxy 
color and the average age of the stellar population. In any case, the RMS of 0.07 dex in R23 will 
often be on the same order as, or even less than the statistical measurement uncertainties on the 
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strong line equivalent widths in high-redshift spectroscopic surveys, even when the signal-to-noise 
of the emission line equivalent widths is as low as 8:1. 

Figure 6 shows a comparison of EWR23 with R^ 3 and R23 for the K92+ and KISS galaxy 
samples. The left column shows the comparison of EWR23 with R^ and the associated residuals 
from the 1-to-l correspondence. There is a good correlation between EWR23 with R^ in the 
upper left panel. The RMS dispersion from the 1-to-l relation is o-(log[i?2 3 ]) = 0.11 dex for the 
combined K92+KISS samples. The left column shows systematic residuals with galaxy color and 
line strength indicating that reddening is significant for some galaxies. The right column again 
shows the comparison of EWR23 with R23 and the associated residuals as a function of equivalent 
width and galaxy color. The residuals for the KISS galaxies are now slightly smaller and much less 
systematic after application of a reddening correction. The RMS dispersion from the 1-to-l relation 
is a (log [-R23]) = 0.09 dex for the combined K92+KISS samples. The residuals are larger than for 
the K92+ and NFGS objects, and the KISS galaxies do not show the systematics with galaxy color 
or line strengths and ratios which the NFGS galaxies exhibit. The lack of systematics with color is 
attributable to slit-vs- integrated light issues, but this does not explain the larger scatter or lack of 
systematics with pure line measures. This latter is possibly a result of the relative homogeneity of 
the KISS sample, as KISS galaxies have stronger emission lines and do no include the diversity of 
more quiescent galaxies with older stellar populations found in the NFGS. 



5.3.1. Effects of Stellar Balmer Absorption 

Ideally, the quantity R23 from which a metallicity is derived should be computed using an 
H(3 line strength which has been corrected for both interstellar reddening and absorption by at- 
mospheres of the underlying stellar population. In practice, the amount of underlying absorption 
is difficult to measure even under ideal circumstances with high signal-to-noise data. Spectra of 
distant galaxies frequently lack the signal-to-noise necessary to measure multiple Balmer lines and 
correct simultaneously for extinction and Balmer absorption in a self-consistent fashion. For galax- 
ies with strong emission lines due to active star formation (i.e., Wjip > 25 A), a correction of a 
few A to the HP line will have a small impact on the derived R23 or EWR23- However, in galax- 
ies dominated by older stellar populations with weak emission lines, R23 or EWR23 will depend 
sensitively on the correction for Balmer absorption. 

Until this point in the analysis, we have not made any corrections for stellar Balmer absorption. 
The effect of underlying Balmer absorption (specifically the amount of absorption in the H/3 line, 
W}{j3{abs) will depress the measured Fjjb an d Wh/3- This leads to systematically large R23 or 
EWR23, an d systematically low oxygen abundances for objects on the upper (metal-rich) branch 
of the empirical calibrations. 

The impact of stellar absorption can be assessed using Figure 7. Using the K92+ galaxy 
spectra, we measured several Balmer lines and performed a self-consistent reddening and stellar 
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absorption solution for each galaxy (see Kobulnicky & Skillman, 1996, for additional details). The 
results indicated a range in underlying H(3 absorption of ~l-3 A with an average value near 2 
A. From these data we calculate a quantity which includes the corrections for reddening and 
stellar absorption calculated for each galaxy. The upper left panel of Figure 7 compares the "raw" 
EWR23 ratio with the quantity R^; the correlation is modest, with a dispersion of o-(-R^ 3 ) = 0.14 
dex and a systematic offset of 0.06 dex. The lower panels illustrate the nature of the residuals as a 
function of galaxy color and line strengths. As might be expected, galaxies with the lowest Wh/3 are 
the most deviant, while galaxies with Wnp > 20 show a much smaller dispersion. A more logical 
approach is also to add a correction to Wjj/3 for the underlying stellar absorption, forming a new 
quantity EWR^- A 2 A correction to Wh/3 was chosen because it was the mean correction needed 
to produce R^ 3 , and is consistent with mean corrections found for other galaxies (e.g., McCall et al. 
1985; Olofsson 1995). The right top panel shows EWR^ 3 plotted against R^. The residuals arc 
now much smaller with cr^R^) = 0.08 dex and a systematic offset of less than 0.01 dex. In the 
absence of direct measurements of the Balmer absorption due to stellar populations, application of 
a 2 A blanket correction to Wh/3 appears to be prudent for the types of starforming galaxies in the 
K92+, DGSS and KISS samples. Samples including galaxies with older stellar populations may 
require larger corrections. 



6. Discussion and Conclusions 



The ratios of the equivalent widths of strong oxygen and hydrogen emission lines from the 
ionized component of distant galaxies can be used as a measure of the global ISM metallicity via 
the substitution of EWR23 fo r -^23- We recommend the use of EWR^ where the Wup has been 
corrected for stellar Balmer absorption assuming a mean correction of 2 A. The typical disper- 
sion from the 1-to-l relation between either EWR23 or EWR^ 3 and the canonical reddening and 
absorption-corrected R23 ratio is cr(i?23) = 0.08 dex. Residuals are somewhat smaller (cr = 0.05 
dex) for galaxies with the largest emission- line equivalent widths (i.e., those having the largest 
rates of star formation per unit luminosity). The translation of these errors into errors in oxygen 
abundance will depend on the value of R23, but for the upper-metallicity branch away from the 
"turnaround" point, the slope of log(0/H) vs log(i?23) is ~l-4 (Pilyugin 2001), so the addition 
error in log(0/H) would be ~0.11 dex. Thus, the additional uncertainty introduced by translating 
a set of measured equivalent widths into the traditional R23 flux ratio diagnostic is comparable to 
or less than the typical observational line measurement uncertainties and systematic errors in the 
R23 to O/H calibration which run 0.15 dex in O/H (e.g., Kobulnicky, Kennicutt & Pizagno 1999 
for a more detailed discussion of the error budget). We anticipate that the method tested here will 
be useful for performing rough chemical abundance estimates in large high-redshift galaxy samples. 
The approach described here will be most useful in a statistical sense when large numbers of objects 
are available for study. Possible applications include understanding the overall chemical evolution 
of star forming galaxies over large intervals of cosmic time (Kobulnicky et al. 2003) or assessing 
the impact of cluster environment and the intracluster medium on the chemical properties of the 
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ISM within galaxies (e.g., Skillman et al. 1996). 
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Fig. 1. — Comparison of dereddened emission line fluxes and flux ratios to emission line equiv- 
alent widths and EW ratios for galaxies from Kennicutt (1992a,b) and Jansen et al. (2000a, b; 
NFGS). Upper panels show the [O IIJ/H/3 flux ratios and EW ratios with a line illustrating the 
1-to-l correspondence. Lower panels show residuals from the 1-to-l relation as a function of equiv- 
alent width and galaxy color. There is generally a strong correlation between flux ratios and 
EW ratios. Panels showing systematic residuals are discussed in the text. Note that the ratio 
(W[oih/Whp)/(I[oii\/IH0) is a measure of l/a. 
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Fig. 2. — Comparison of dereddened emission line fluxes and flux ratios to emission line equivalent 
widths and EW ratios for galaxies from the K92+ galaxies and the KISS (Salzer et al. 2001) 
emission line galaxy survey. Upper panels show the [O IIJ/H/3 flux ratios and EW ratios with the 
line illustrating the 1-to-l correspondence. Lower panels show residuals from the 1-to-l relation as 
a function of equivalent width and galaxy color. There is generally a strong correlation between 
flux ratios and EW ratios. Panels showing systematic residuals are discussed in the text. 
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Fig. 3. — Comparison of the ionization parameter quantity W\p iu]/W\o il\ with F\q jjn /F\q jn 
and the dereddened ratio I\o iin/ho II\ f° r the K92+ and NFGS galaxies. Lower panels show 
the residuals from the 1-to-l correspondence as a function of line strength and galaxy color. The 
tight correlation between (W[q iii]/W[o ii])/(F\o ni]/F[0 an d (B — V) is expected, as the latter 
simply represents the ratio of spectroscopically measured continuum flux under the lines. The ratio 
(W[oiii]/W [OII] )/{I [OIII] /I [OII] ) gives the unknown value a. 
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Fig. 4. — Same as Figure 3 except with KISS in place of NFGS galaxies. The lack of a tight 
correlation between (B — V) and (W\o ni]/W[o in) /(Ftp nn/ 'Ftp in) is most likely due to observing 
different regions of each galaxy for spectroscopy and photometry (see text). 
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Fig. 5. — Comparison of the quantity R23 and with EWR23- R23 * s -^23 without correction 
for reddening. The strong correlation between R^ 3 and EWR23 (upper left panel) is even stronger 
for R23 and EWR23, suggesting that oxygen abundances can be estimated from equivalent width 
ratios as well as from dereddened line fluxes. The improvement of the correlation with R23 shows 
that EWR23 is successful in correcting at least part of the reddening. The RMS dispersion from the 
1-to-l relation is o"(log[i?23]) = 0.07 dex. Lower panels explore the residuals in the correlation as a 
function of line strength, line ratio, and galaxy color, left column: no correction for extinction or 
underlying Balmer absorption in either quantity; right column: line fluxes have been corrected for 
reddening. There remains a small systematic trend with color, but, excluding only a few outliers, 
the EWR23 an d R23 values agree to about ±0.1 dex for this sample. 
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Fig. 6. — Comparison of the quantity R23 and with EWR23 f° r the K92+ and KISS galaxy 
samples. Upper panels show the excellent correlation between R23 and R^ with EWR23- Lower 
panels show the residuals from the 1-to-l correspondence as a function of line strength and galaxy 
color. Residuals are larger, but less systematic, for the KISS galaxies than for the NFGS and K92+ 
galaxies. 
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Fig. 7. — Effects of correcting for underlying H(3 absorption in K92+ galaxies. The quantities 
#23 and EWR2Z are analogous to -R23 and EWR23 but with the H(3 measurements corrected for 
stellar absorption by 2 A, an average value found in a self-consistent fit to both the reddening and 
stellar absorption using the available Balmer lines (see text). Upper panels show the correlation 
between -R23+ and with EWR23 and EWR 2 s+- Lower panels show the residuals from the 1-to-l 
correspondence as a function of line strength and galaxy color. Galaxies with Wh/3 < 15 A are 
most seriously affected by the lack of correction for stellar absorption. 



